Quantum effect of inductance on geometric Cooper-pair transport 



(N 

o 

(N 

> 
O 

o 

(N 



43" 

-I— > ■ 



> 

in 
in 

en; 

ON' 

o: 



X 



Juha Salmilehto 1 and Mikko Mottonen 1, 2 

1 COMP Centre of Excellence, Department of Applied Physics, 
Aalto University, P.O. Box 13500, FI-00016 AALTO, Finland 
2 Low Temperature Laboratory, Aalto University, P.O. Box 13500, FI-00076 AALTO, Finland 

We introduce a model for a flux-assisted Cooper-pair pump, the sluice, which is used to study 
geometric charge transport. Our model allows for a nonvanishing loop inductance going beyond the 
usual treatment with an exact phase bias. We derive the device Hamiltonian and current operators 
for different elements of the system and calculate the pumped charge carried by the ground state in 
the adiabatic limit. We show that extending the model beyond the exact phase bias, has a weak but 
potentially non-negligible effect on the charge transport. This effect is observed to depend on the 
external flux bias. The adiabatic energy level and eigenstate structures are studied and the unusual 
features are explained. Finally, we discuss the requirements of adiabaticity. 
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I. INTRODUCTION 



The methodology of constructing superconducting 
Josephson junction circuits^ has provided a wondrous 
test-bed for many quantum phenomena. Perhaps one 
of the most intriguing applications is the possibility to 
study physics related to quantum phases 2 ., most no- 
tably the Berry phase 3 -! 4 -. Understanding the properties 
of these phases 5 - - — has given rise to geometric quantum 
computin g 10 ! 11 where the robustness against a certain 
type of noise is an inherent property of the computing 
scheme. 

One of the most extensively studied Josephson devices 
related to geometric phases is the Cooper-pair sluic o 12 ' 13 . 
It builds on the idea of implementing an array of Joseph- 
son junctions where single Cooper pairs are adiabati- 
cally transported in a controlled manner— ~— establish- 
ing a link between the pumped charge and the Berry 
phas o 18 ! 19 . The introduction of superconducting quan- 
tum interference devices (SQUIDs) operating as tun- 
able junctions adds to the pump control^ 2 -. Cooper- 
pair pumps based on such array structures have been 
theoretically shown to produce very accurate charge 
quantization 2 ^— as well as having the potential for im- 
plementing non-Abelian structure a 23 ' 24 . For the sluice, 
the connection between the pumped charge and the 
Berry phase has been shown both theoreticall y 18 ! 25 and 
experimentally 2 ^. Furthermore, the development of the 
sluice has given rise to proposals to implement geometric 
quantum computing featuring similar architectures 2 ^ 2 ^. 
The sluice itself has been proposed to be used as a 
Landau-Zener-Stiickelberg interferometer—. 

During the past few years, theoretical effort has been 
put into modeling the sluice coupled to a dissipative 
environment. Using the usual reduced-density-operator 
method a 30 i 31 , such a coupling leads to charge noncon- 
servation. Recent papers have dealt with this problem 
using either superadiabatic basest— or the Floquet 
theor y 36 i 37 . Both methods also have the potential for 
describing pumping when the driving is nonadiabatic, a 



regime recently studied experimentally 3 ^,. Furthermore, 
a current induced directly by the environment has been 
shown to emerge for certain types of noise operators 3 ^ 
leading to a development of a conservation law for all 
operator currents in open quantum systems^ .. To date, 
all considerations of the sluice have assumed an exact 
phase bias. Namely, the inductance of a superconduct- 
ing loop at the heart of the device has been assumed 
so small that the total gauge-invariant phase difference 
across all the weak links along the loop is fixed exactly 
by the penetrating magnetix flux. This assumption has 
simplified the theoretical analysis which has still been 
shown to provide predictions in agreement with the im- 
plemented experiments 2 ^. However, these experiments 
have not been carried out with very high precision and 
given that the sluice has been widely discussed as a can- 
didate for a metrological current source, it is of utmost 
importance to study whether the finite realistic loop in- 
ductance plays a role here. 

In this paper, we introduce a theoretical model for the 
sluice which allows the inclusion of a nonvanishing loop 
inductance. The derived Hamiltonian accounts for the 
emerging additional quantum degrees of freedom. The 
current operators for the superconducting island and the 
effective inductor are derived from first principles. We 
calculate the geometric charge carried by the adiabatic 
ground state emphasizing the features due to the finite 
loop inductance and show that the results given by the 
model with vanishing inductance are valid to a rather 
good degree with experimentally relevant parameters. 
Variation of the additional parameters included in our 
model is shown to weakly affect the pumped charge. We 
study the instantaneous eigenenergies and eigenstates to 
point out the differences between our more general model 
and the case of the exact phase bias. Finally, we consider 
the adiabaticity criteria for the system. 

The paper is organized as follows. In Sec. Ill Al we 
introduce our model for the sluice and derive the Hamil- 
tonian for it. In Sec. Ill Bl we derive current operators 
for the different parts of the system without resorting to 
operator derivatives. In Sec. IIII Al we study analytically 
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the transferred geometric and dynamic charges in first or- 
der perturbation theory. The charges are found equal to 
the ones obtained using the original model. In Sec. MI Bl 
we perform a numerical analysis of both the system and 
its charge transport features in the adiabatic limit. The 
geometric charge carried by the ground state is shown to 
differ slightly from the original model. The features of 
the energy level structure and eigenstates are explained 
and the requirements of adiabaticity are discussed briefly. 
We conclude the paper in Sec. IIVI 

II. MODEL DEVICE 
A. Hamiltonian 

The device structure we use is an adaptation of the 
original model used for the Cooper pair sluice^. The 
original model consists of a superconducting island iso- 
lated by two SQUIDs positioned in a large superconduct- 
ing loop [see Fig. [1] for our model resembling the origi- 
nal structure with respect to the device operation]. The 
island is capacitively coupled to a gate terminal used 
to manipulate its charge state. The leads are taken to 
have capacitive couplings to a ground terminal. The 
self-inductances of the SQUID loops are assumed neg- 
ligible so that they act as tunable Josephson junctions 
whose Josephson energies can be controlled with the ex- 
ternal magnetic fluxes threading the loops. Phase-biasing 
is typically introduced for the device by assuming that 
the total inductance of the large loop is negligible imply- 
ing that the total superconducting phase difference over 
all the weak links in the loop is a real number exactly 
determined by the penetrating magnetic flux. This as- 
sumption allows one to determine the system dynamics 
using only the island charge basis. 

We improve on the existing model by introducing an in- 
ductive element allowing us to depict, among other induc- 
tance sources, the self-inductance of the loorjii. Ground 
couplings are established from both sides of the element 
using two identical shunt capacitors. The resulting de- 
vice is shown in Fig. [TJ The Hamiltonian for the device 
can be written as [see Appendix [S] for details] 

E 

H = E c (n~ n g f + E v fi 2 + H Fm + -^(<p- f> ) 2 . . 
- E Jr cos(0 + ip/ 2) - Eji cos(0 - ip/2), 

where (f> is the island phase operator, tp is the opera- 
tor for the total phase difference over the SQUIDs, n is 
the Coope- pair number (CPN) operator for the island 
and n = (ni — n r )/2 is the Cooper-pair number opera- 
tor for the feed; that is, it describes the average charge 
imbalance between the left and right leads. The CPN 
operators for the left and right leads are given by ni and 
n r , respectively. To allow for a form resembling the case 
of the exact phase bias, we omitted explicitly writing all 
identity operators but will reestablish their use in the 




FIG. 1. Circuit diagram of the Cooper pair sluice. The exter- 
nal fluxes threading the left and right SQUIDs are denoted by 
and <E> r , respectively, and $ denotes the operator for the 
total flux threading the large superconducting loop. The op- 
erators for the phase differences over the SQUIDs are marked 
by ipi and (p r and they are defined in the direction specified 
by the arrows in the figure. The internal capacitances of the 
SQUIDs are Ci and CV, the loop is connected to the ground 
terminal via two identical capacitors with capacitances C v , 
and L — Lq + Lk denotes the total inductance of the loop 
comprising of both the geometric Lq and the kinetic Lk con- 
tributions. The gate capacitance C g is used to manipulate the 
island charge with the gate voltage V g . The operators h, hi, 
and h r are for the Cooper-pair numbers of the corresponding 
islands. Note that the superconducting islands described by 
hi and h r are inductively coupled. The choice of the charge 
orientation at the capacitors for positive voltages across the 
capacitor plates is marked by {+, — } including the effective 
junction capacitors denoted by dashed lines. 



following presentations of the Hamiltonian. Here, Eq is 
the charging energy scale for the island [see Eq. (|A2|) ]. E v 
is the charging energy scale for the feed [see Eq. (|A2[) ]. 

El = j- (f 2 -) is the inductive energy scale due to the 
finite total inductance L of the loop, Hpm is an addi- 
tional charging Hamiltonian stemming from the deriva- 
tion [see Eqs. (TM} and ([A4])]. <p = 2 ^ axt is related the 

external magnetic flux $ C xt, % = 9 2e 9 is the normalized 
gate charge, and Eji,Eji are the Josephson energies of 
the tunable junctions. The flux quantum is denoted by 
$0 ~ 2.07 fWb, the elementary charge is denoted by e 
and L = Lq + Lk , where Lq is the total geometric in- 
ductance of the loop and Lk is the corresponding kinetic 
inductance. The derivation of the Hamiltonian and the 
exact definitions of all the quantities are given in Ap- 
pendix IaI 

We have written the Hamiltonian using two pairs of 
conjugate variables; that is, (n, c/>) and (n, (p). Especially, 
the definition of the feed CPN operator allows us to treat 
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the inductive part of the Hamiltonian. We identify 
*2 , El 



H 



v - E v n + —(cp - ip ) 

d 2 E L 2 



(2) 



as being the Hamiltonian for a phase-shifted quantum 
harmonic oscillator, for which the effective mass is m = 
h 2 /(2E v ) and the effective angular frequency is u) = 
v /2E v E L /h. If we denote d = {2E v /E L ) 1 / i , the cor- 
responding eigenproblem H v \^>k) = E^ \^k) defines the 
Fock states for the feed 
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x H k 



<£ ~ Po 
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(3) 



where k is any non-negative integer and Hk is the fcth 
Hermite polynomial. The corresponding energies are 
given by 



E k = ^2E~E~ L (k + 



(4) 



Notice that the canonical commutation relation is 
[<p, n] = i implying that the effective momentum oper- 
ator is hh. Defining the relevant bosonic creation and 
annihilation operators for the feed Fock states as a' and 
a, respectively, allows us to write the total phase differ- 
ence operator as 



<P= \ ~2E~ J (a + a'j + ^o- 



(5) 



Using the above-mentioned notation, the Hamiltonian 
can be written as 



placement. This notation is typically utilized in quantum 
optics. Finally, we rewrite the Hamiltonian as 



H = H, 



FIN 



(8) 



+ H B ® [D(igo) - %>] + H B ® [D\ig ) - % 

where H$ = Ec(h—n g ) 2 —Ej r cos(4>+ipo/2)—Eji cos(0— 
tpo/2) is the Hamiltonian for the original model with 

exact phase bias and Hb = —^-^ — [Ej r h+ + Ejih-]. 
Note that if the total inductance is negligible (L — >• 0, 
El — > oo), the system dynamics are accurately described 
by the original Hamiltonian acting only on the island de- 
grees of freedom. 



B. System currents 

To study the currents flowing in the device, we describe 
the temporal change of the expectation value of the feed 
charge operator using the corresponding current operator 
I = — ^jf[H, h], where we exploited the usual current 
identity applicable for any arbitrary charge operator that 
does not depend on time explicitl y 25 ! 35 . We thus proceed 
to write the feed CPN operator as 



n = —i (a — ah, 

45o 



(9) 



which allows us to write the general form for the current 
operator as 



2hg 



[if, I„<8> (a -a*)]. 



(10) 



Exploiting the Hamiltonian derived in the previous sec- 
tion and utilizing some of the properties of the dis- 
placement operator familiar from the coherent state 
formalism^, the current operator becomes 



H = H n ® I v + I n <g> H v + H Fm 



Ej r 



Ej 



{h+ ® e lv °/ 2 D(ig ) + h- <g> e - tVo/2 D(-ig )} 



-^{n+ ® e-^°l 2 b{~ig Q ) + n_ <g> e^°/ 2 D{ig )} , 



(6) 



where H„ = Ec(n — n g ) 2 , I„ is the identity operator in 
the island CPN basis, I„ is the identity operator in the 
feed Fock basis, h + = exp[i0] (n_ = exp[— i<p]) is the 
raising (lowering) operator in the island charge basis and 
D(igo) — cxp[igo(a + a')], where 



( e v y/ 4 

50 \V2E L ) 



(7) 



The operator D(ig a ) is the displacement operator for one 
mode in the quantum phase space with igo being the dis- 



I=j[H B ® D(igo) -H B ® D\igo)] 



2e 
~W 



E V E\ 



1/4 



I„(g) (a + a f ). 



(11) 



Note that the first and third terms of the Hamiltonian in 
Eq. (JSJ) do not contribute anything to the current opera- 
tor due to commutation. If the total inductance is neg- 
ligible, the first term on the right-hand side in Eq. (TTTj) 
yields the current operator used in the original model. 

The current operator in Eq. pip can be interpreted in 
terms of currents flowing in the system. The operator / 
describes the temporal change in the charge imbalance 
and, hence, the transfer of charge between the shunt ca- 
pacitors or the effective junction capacitors. Due to the 
device structure, the charge can be transferred via two 
separate routes: through the inductor or through the is- 
land, and the total transfer must be described by /. 

To quantify the argument given above, let us present 
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FIG. 2. Schematic diagram of the electric currents for the 
model used to analyze the Cooper-pair sluice. See text for 
further explanation of the current operators. 



the current schematic for the sluice in Fig. [2] In the fig- 
ure, I j is the operator for the average current tunneling 
through the junctions, II is the operator for the cur- 
rent flowing through the inductor and 1\ (I2) describes 
the current flowing to the junction system from the left 
(right). Furthermore, we denote the operator for the cur- 
rent leaving the left shunt capacitor by and for the cur- 
rent arriving to the right shunt capacitor by I r . By the 
definition of the feed charge operator presented in Ap- 
pendix [Aj we have I = Ii+Ij — Ii, where we use the fact 
that the system of islands is capacitively isolated from 
any external current sources. By current conservation, 
ii = J; + II and, hence, I = II + Ij. Thus, the simple 
argument given above about the separation of / into dif- 
ferent current operators for the system has been shown 
valid. In addition, we have —I r = II + Ii- If the current 
arriving to the island equals the average current tunnel- 
ing through the junctions I2 = lj, then the conservation 
of total charge = implies that I = —I r = —Ii in- 
dicating that the current operator describes discharging 
(charging) of the right (left) shunt capacitor. 

Identifying the current operators for the inductive el- 
ement and for the tunneling current across the island is 
straightforward. As described in Appendix A, the induc- 
tive current can be derived by studying the operator for 
the total flux penetrating the loop $ = $ cxt + LqIl and 
its relation to the operator for the total phase difference 



over the SQUIDs 



*0 
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inductive current operator given in Eq. 



Rewriting the 
fA5]) yields 



II 



2e , 



<Po) 



2e 8H im 
h dip 



(12) 



where the last form is just to show that the usual manner 
of defining via operator derivates is also applicable. For- 
mulating this current operator in terms of the creation 



and annihilation operators allows us to write 
2eE Ll 



' : ^Y 1 [2g li (a + a)) + ip - <p ] 
(27 i «®(« + « t )' 



2e 



(13) 



which can be identified as the second term in I in 
Eq. (jlip . Using the operator derivative, we could also 
identify the first term in Eq. as stemming from the 
Josephson part of the Hamiltonian and, hence, describ- 
ing the charge passing through the island. However, due 
to the above-derived separation of currents, such appli- 
cation of operator derivatives is not necessary in order to 
define the island current operator. 

In the adiabatic limit, there is no voltage over the in- 
ductor so that the potentials on the left and right islands 
are equal. Thus, there is no potential difference that 
would induce a change in the charge imbalance between 
the left and right leads, and hence (/) = at all times; 
that is, there is no net current flowing in or out of the 
shunt capacitors. In this case, we can define the charge 
circulating in the loop to be described by either II or Ij . 



III. ANALYSIS 

Using the assumption of negligible loop inductance, the 
geometric charge transferred adiabatically through the 
sluice during one parameter cycle has been theoretically 
studied in several reference s 12 ' 25 ' 32 ' 33 ' 35-37 ' 39 . However, a 
model enabling one to use a non-negligible inductance for 
such transfer has not been studied previously. We write 
the geometric charge transferred by the mth adiabatic 
state in the usual manne r 12 ! 15 ' 25 



n (G,m) 



-2fSsm 



(m(q)\i {I/L) (q)\k(q)) 



E 



x (k{q)\V f\m(q)) 



n{q) 
dq, 



E k (q) 



(14) 



where \i(q)) is the ith instantaneous eigenstate of the 
Hamiltonian in Eq. ([5]) and Ei(q) is the corresponding 
eigenenergy. We present everything as a function of the 
point in the control parameter space denoted by the vec- 
tor q whose scalar components are the external control 
parameters. The integral is over a contour 7 in the pa- 
rameter space. In the case of cyclic evolution, the con- 
tour is closed and the pumped charge is related to the 
Berry phase obtained during the time evolutio n 18 ' 25 ' 26 . 
The pumped charge is fully determined by the contour 
and, hence, it does not depend on the pumping speed. 
For analytical calculations, we also define the dynamic 
charge carried by the supercurrent as 



Q 



(D,m) 
(I/L) 



1// 



dt(m(t)\i (I/L) (t)\m(t)), (15) 
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where the pumping frequency is denoted by /. 

In the following, we assume that the island is capac- 
itively symmetric Ci = C r = Cj and that C g <C C v . 
According to Appendix |A1 //fin can consequently be ne- 
glected in the total Hamiltonian. Furthermore, we as- 
sume that the external control parameters are manipu- 
lated slowly enough for the evolution to remain adiabatic. 



A. Perturbative approach in the lowest order in go 

Perturbation theory can be applied in the lowest order 
in 90 to obtain analytical information on the properties 
of the system. Decomposing the displacement operator 
yields 

D(ig ) =l v +ig a {a + a))+0{gl). (16) 

If E<p <C El, go is small and the linear term suffices so 
that the Hamiltonian in Eq. J5| is accurately approxi- 
mated by 

H ss H s ® l v + l n ® H v + ig (H B - H B ) ® (a + a f ). 

(17) 

Similarly, taking only the linear term in Eq. (|lip yields 
7 = I s ® % - g ^{H B + H B ) ®{a + a)) 

+ -H-y-M I«® (a + at), 

where we have identified Is = ^(H B —H B ) as the current 
operator in the original model. For future reference, we 
denote the first term on the right-hand side as the zeroth 
order island current operator jj ^ and the second term as 
the first order island current operator /j 1 ^ . 

After an analytical calculation, it can be shown that 
no first-order perturbative correction in go emerges to 
the dynamic or geometric currents across the island com- 
pared with the original model. Furthermore, the corre- 
sponding currents through the inductor can be shown to 
be equal but of opposite sign up to this order. This obser- 
vation agrees with our previous assumption that the total 
temporal change of the expectation value of the feed CPN 
operator is zero in the adiabatic limit. See Appendix 151 
for details of the perturbative calculation. Exploiting 
higher-order perturbation theory would require a higher- 
order expansion of the displacement operator in Eq. (|16[) 
leading to an increasingly more difficult analytical anal- 
ysis. Furthermore, we wish to study a region where go is 
significantly large so that the perturbative approach is no 
longer practical. Hence, we resort to numerical analysis 
below to study the current contribution arising from the 
finite loop inductance. 



B. Numerical analysis 

We begin our analysis by discussing the physical en- 
ergy scales. Assuming that C g is nearly negligible, the 
charging energy scales are determined by Eq ~ e 2 /Cj 
and E v w 4e 2 /(C ¥ , -I- Cj) [see Appendix |A"]. To establish 
a comparison with recent studie o 32 ' 33 > 35 ~ 37 i 39 , we assume 
that Ec/ks = 1 K, where k B is the Boltzmann constant. 
This defines the junction capacitances Cj = Ci = C r to 
be of the order of femtofarads. We study a range from 
0.05 to 500 fF in the shunt capacitances C v . We can as- 
sume, for example, that the shunt capacitances are pri- 
marily due to capacitive couplings to a bottom ground 
terminal via bonding pads. If the device is constructed 
on a 500 /im thick silicon wafer with square-shaped pads, 
the equivalent range in the length of the side of the pads 
is from 15 /um to 1.5 mm yielding experimentally rel- 
evant values. Using the above-mentioned capacitances, 
the maximum feed energy scale is limited by Cj so that 
Etp decreases from 4Ec to 10~ 2 Ec with increasing pad 
capacitance. 

The inductive energy scale is determined by El = 

X (§ tt ) [ see Appendix IA"] . We approximate the geomet- 
ric inductance of the wire by Lq ps jil. where [i is the 
permeability of the wire material and I is the length 
of the large superconducting loop. The kinetic induc- 
tance can be written in the BCS formalism as Ik = 

(™) S tanh(A/2fc fl T) > wliere W iS the width ° f tlie SU P er " 

conducting wire, R sq is the sheet resistance of the wire 
in the non-superconducting state, A is the superconduc- 
tor energy gap and T is the wire temperatur o 43 i 44 . We 
assume that TcT c , where T c is the critical tempera- 
ture of the superconductor, and that the wire is made of 
aluminum such that w is of the order of a micrometer—. 
In addition, we assume that the thickness of the wire is 
of the order of tens of nanometers or larger limiting the 
kinetic inductance to be of the same order than the geo- 
metric one or smaller. The feasible length scales for the 
loop range between 10 ^m and 1 cm corresponding to 
values between l0 3 Ec and Eq in the inductive energy 
scale. 

The standard control parameter cycle utilized in 
Refs. HU, HH, I35l - l37l . and 113 is adopted and the device is 
operated in the charging regime Ec 3> max{ Eji, Ej r }. 
The parameter cycle is presented in Fig. [3] [see Ap- 
pendix [C] for a detailed description]. Furthermore, the 
normalized gate charge n g is manipulated near 1/2 so 
that the required island states are limited near this value. 
In the original model, using only the two closest island 
CPN states \n = 0) and \n = 1) was sufficient. We con- 
duct our simulations using a basis with an even num- 
ber Ni of island CPN states such that the basis con- 
sists of Ni/2 pairs of states \n = 1 — i) and \n = i), where 
i = 1, . . . , Ni/2, and Nf feed Fock states with the lowest 
energies. 

According to Eq. (JSJ , the interesting features are char- 
acterized by two factors: the energy gap Af m ^/IE^El 
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FIG. 3. (Color online) Time dependence of the control pa- 
rameters Eji, Ej r , and n g during a Cooper-pair pumping 
cycle. 



between the states in different Fock subspaces and go = 
{ e \ 1/4 

( 32 e l ) dictating the strength of the coupling between 
the island and feed degrees of freedom. To assess the 
implications of these features with respect to the original 
model, we consider three different points in the (El, E v )- 
space. Firstly, we establish a base point corresponding 
roughly to the energy scales used in the experiments of 
Ref. HI using E L = 10 2 E C and E v = E c . Secondly, we 
depict the point of smallest Af within the feasible en- 
ergy scales using El = Eq and E v = 10~ 2 Ec- Thirdly, 
we illustrate the point of largest go using El = Ec and 
E v = AEc- We perform the simulations for all values of 
the flux bias and retrieve a pumped charge close to the 
usual cosine dependence Q G = 2e(l — 2ecos<y9 ) + 0(e 2 ), 
where e = EJ^/EJg*, k e {l,r}^.. Identical SQUIDs 
are assumed and we exploit the parameter cycle given in 
Appendix [C] The flux bias dependence of the pumped 
charge at the three selected parameter points is given 
in Fig. HJa) and the charge differences between the base 
point and the other two points in Fig. HJb). The charge 
transferred at the base point is equal to the one given by 
the original model within the accuracy of the simulations. 
Note however that the accuracy of both simulations is or- 
ders of magnitude lower than the metrological accuracy 
10~ 8 . Furthermore, the charges transferred through the 
island and through the inductor were found equal in all 
the simulations validating the assumptions made about 
the adiabatic limit in Sec. Ill Bl 

Near the phase values Lpo = tt/2 + Ntt, where N € Z, 
the pumped charge is robust against variation of the ad- 
ditional parameters included in our model whereas at 
(fio = Ntt, the absolute difference between the charges 
reaches its local maxima; see Fig. [31(b) . Furthermore, 
averaging over the flux bias yields a vanishing charge 
difference. Both of these features are likely caused by 



the shape and symmetry of the parameter cycle. The 
two quantum degrees of freedom in our model are cou- 
pled by the tunneling term in the Hamiltonian. A crude 
estimate would thus state that the feed degree of free- 
dom acts as if the original model obtains a reshaping 
of the parameter cycle for the SQUIDs explaining the 
vanishing effect on the pumped charge at approximately 
tpo = f/2 + Ntt, the same point where the vanishing ef- 
fect of the residual Josephson energies appeared in the 
original analytical model. Note that treating the induc- 
tance classically within the original model fails to capture 
the above-mentioned dependence on the external flux^. 
The observed weak influence of altering the parameters is 
especially remarkable since, according to Eq. (fT4")l . the ge- 
ometric charge is not affected only by changes in the state 
carrying it but also by changes in all the other states. 

We study the instantaneous energy level structure and 
eigenstates emerging from our model focusing near the 
point of maximum robustness; that is, we select the ex- 
ternal magnetic field to yield tpo = tt/2. According to 
the original model, this point is of special interest as the 
lowest-order effect of the residual Josephson energy on 
the pumped charge vanishes. The energy level structure 
at the different parameter points is given in Figs.[5fa-c). 
For the probability distribution of the ground-state pop- 
ulation to different basis states P(n, k), where n refers to 
the island CPN state and k to the feed Fock state, we only 
show the elements that exceed 10 -2 at any time instant 
in Fig.[3Jd). At the base parameter point, there is hardly 
any mixing between different feed Fock subspaces which 
are separated by a large energy gap. Restricted to one 
such subspace, the populations resemble the ones given 
by the original model. At the second parameter point, 
the subspaces are no longer separated by a large gap and 
the energy diagram is greatly affected. Especially, the 
higher energy levels obtain a complicated structure hav- 
ing multiple avoided crossings where the energy gap is 
very small during the cycle. The ground-state distri- 
bution is still similar to that in the previous case, but 
the other adiabatic states change significantly. We give 
the first excited state in Fig. [5] showing that the state 
remains isolated to the first excited feed Fock subspace 
during the time evolution except in the vicinity of the 
avoided crossings where some of the population leaks to 
the other feed Fock subspaces. Note especially that the 
first excited state is no longer confined to the same Fock 
subspace as the ground state which may be significant 
for nonadiabatic and dissipative dynamics. At the third 
parameter point, the parameter go is more than 4 times 
that at the base point. However, the effect of this increase 
on the level structure seems small and the most signif- 
icant changes, including the emergence of the avoided 
crossings at the higher levels, can be explained by the di- 
minishing of Af . The ground-state distributions for the 
third parameter point and the base point are very simi- 
lar [see Fig. EJd) for the distribution at the base point], 
but a closer inspection reveals a small deviation between 
the two near the degeneracy points. This was to be ex- 
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FIG. 4. (Color online) (a) Pumped charges and (b) charge differences as a function of the scaled external magnetic flux bias 
fo = 22f^t. We use E c /k B = 1 K, EJf/Ec = EJT/E a = 0.1, Ej; n /E c = Efl n /E D = 0.003, < ax = 0.8, nf n = 0.2. The 
parameter points in (a) are El = 10 2 Ec and E v = Ec (base point), El = Ec and E v — W~ 2 Ec (point 2), and El = Ec and 
E v — 4Ec (point 3). The dashed line depicts the approximate analytical result derived for the original model and the dashed 
rectangle indicates the area of the inset. In (b), we give the geometric charge differences between the base point and point 2, 
and between the base point and point 3 (from bottom to top at <^>o = tt). For the simulations, we use Ni — 2 and Nj = 25. 



( b ) 1.. 




FIG. 5. (Color online) Energy level structures and ground-state probability distributions of the sluice as a function of scaled 
time tf, where t is time and / is the pumping frequency. Lowest instantaneous energy levels for (a) El = 10 Ec and E v = Ec, 
(b) E L = E c and E v = 10~ 2 E C , and (c) E L = E c and E v = 4E C during a pumping cycle detailed in Appendix [C] The 
penetrating magnetic flux is selected so that <^o = w/2 and the other parameters are equal to the ones used in Fig. [4] The 
distribution of the ground-state populations P(n,k) on different basis states [see text for further details] is given in panel (d). 
The differences in the distributions between the different parameter points are smaller than the line width in the figure and, 
hence, we only present the distribution for the base point for clarity. We use Ni = 4 and Nf = 25. 
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FIG. 6. (Color online) Distribution of the first-excited-state populations P(n, k) on different basis states [see text for further 
details] for El = Ec and E v = 10 _2 £c as a function of scaled time tf, where t is time and / is the pumping frequency. Other 
physical parameters are the same as in Fig. [4] with tpo = n/2. For the simulations, we use Ni = 4 and Nf = 25. 



pected as the terms involving the displacement operator 
in Eq. © couple to the Josephson terms of the tunable 
junctions. 

Since the robustness of the geometric charge carried 
by the ground state has been shown above, we turn our 
attention to the requirement of energy separation. Both 
Landau-Zener transition o 29 i 38 i 46 and nois o 32 i 35 i 37 i 47 may 
excite the system causing the transferred charge to drop. 
In the original model, the parameter cycle was designed 
such that the two lowest states were well separated from 
the rest of the eigenspectrum. In Fig.[SJa), this assump- 
tion clearly applies for our model, too. However, the 
assumption is broken if Af is of the order of the lowest 
energy gap within a single Fock subspace as shown in 
Fig. [nib) . Hence, using the same driving speed as in the 
original model does not necessarily mean that the transi- 
tions are localized at some well-defined avoided crossings 
similar to the ones in Fig. [5la). Furthermore, the en- 
ergy separation between the higher states is very small 
at the emerging avoided crossings indicating that a com- 
plex process of excitations involving multiple levels may 
occur. This is potentially detrimental for pumping. 



IV. CONCLUSIONS 

We introduced a model for the Cooper pair sluice uti- 
lized to observe the features of geometric charge trans- 
port. Our model allows the inclusion of a nonvanishing 
loop inductance extending the theory beyond the usu- 
ally applied exact-phase-bias picture. The Hamiltonian 
for the device was derived by identifying relevant pairs 
of conjugate operators, and the current operators for the 
different elements of the system were formulated. 

The methodology used to derive the model can be 
straightforwardly adapted to more complicated systems. 
A direct application of the methodology includes the pro- 
posal for geometric quantum computing using a loop with 



two islands separated by SQUIDs^ 7 -. For this device, the 
practical operation window in the extended parameter 
space should be carefully surveyed since the ground state 
degeneracy might be lifted invoking error for each holo- 
nomic gate. In addition, we generally expect that using 
our methodology for modeling dissipative systems yields 
a more exact description of the effect of the environ- 
ment. Applying the full quantum description is especially 
important for systems where the flux noise reveals im- 
portant features, such as the dissipative currents flowing 
in the sluice 3 ^, or where it becomes the dominant noise 
source due to inbuilt robustness against charge noise, as 
is the case for the transmon qubit 4 ^. 

We analyzed the charges carried by the adiabatic states 
during a pumping cycle both analytically and numeri- 
cally. We showed that the difference in the transferred 
charges between our model and the exact-phase-bias pic- 
ture is of the second order or higher in the strength of 
the coupling between the island and feed degrees of free- 
doms. Within experimentally feasible energy scales, we 
retrieve the known cosine dependence of the geometric 
charge with respect to the external magnetic flux with 
good accuracy. The geometrically transferred charge dif- 
fers slightly from the original model when the additional 
quantum degrees of freedom are included using experi- 
mentally relevant values. The effect of altering the addi- 
tional parameters included in our model on the charge de- 
pends on the external magnetic flux bias such that robust 
bias points emerge. The dependence differs from what 
the classical treatment of the loop inductance within the 
original model predicts and, hence, emphasizes the ne- 
cessity of the quantum approach. In addition, applying 
the full quantum picture does not necessarily lead to a 
decreased geometric charge but the charge can also be 
increased by an appropiate selection of the loop induc- 
tance, the shunt capacitances, and the magnetic field. 
Even though the charge is only weakly influenced by the 
changes in the loop inductance, future experiments may 
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be able to probe the variation using experimental setups 
based on current state-of-the-art method a 25 i 26 ' 38 . 

Based on the geometric nature of the pumped 
curren t 21 i 25 , the sluice and its generalizations have some 
potential to be used as quantized current sources for 
metrological applications. However, the metrological 
source requires extreme accuracy such that the maximum 
relative error is ~10~ 8 . The analysis presented in this 
paper shows that applying our model is necessary in the 
studies of the operation and feasibility of the sluice as 
such a source^ 2 -. This is because experimentally relevant 
variations in the additional parameters included in our 
model, especially the self-inductance of the loop, have 
an effect on the pumped charge that is orders of magni- 
tude larger than the relative error required to realize a 
metrological current standard. 

Finally, we studied the instantaneous energy level 
structure and the probability distributions of the instan- 
taneous eigenstates to different basis states near a max- 
imally robust bias point. The energy separation of the 
feed Fock subspaces largely determined the instantaneous 
energies such that a complicated array of avoided cross- 
ings emerged when the separation was small. The ground 
state remained significantly robust whereas the higher 
states were prone to changes. Changing the additional 
energy scales emerging in our model was shown to be 
able to cause significant changes in the system energies. 
In particular, the energy separation between the ground 
state and the first excited state may be diminished and 
temporally altered implying that the time-local Landau- 
Zener transition probability is generally increased and 
the transitions are not necessarily localized at avoided 
crossings similarly to the original model. Drastic adjust- 
ment in the driving speed may be required to experimen- 
tally remain in the adiabatic limit. In addition, reduction 
in the energy separation between the ground state and 
the first excited state potentially causes the environment- 
induced transitions to become more prominent between 
these states. 
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Appendix A: Constructing the device Hamiltonian 

We denote the Cooper pair number (CPN) operators 
for the charge in the tunable SQUID junctions by hi,h r 
and for the charge in the shunt capacitors by n v> i,h v>ri 
where the subindices I and r denote the left and right 
capacitive elements, respectively. The CPN operators 
for the superdonducting islands can then be written as 
n = h r —ni, hi = ni—n^j and h r = —h r +h Vtr as detailed 
in Fig. [TJ In the following, we refer to h as the island 
CPN operator and define n = = - "*■■'+"*■■*■ 

as the feed CPN operator describing the average charge 
imbalance between the superconducting islands that con- 
stitute the left and right leads. The only way to affect 
the imbalance is by charging the shunt capacitor system 
corresponding to charge transport from one capacitor to 
the other or by charging the capacitive junction system 
corresponding to charge transport between the effective 
junction capacitors. 

The system of superconducting islands is electrically 
isolated so that the net charge is constant, i.e., h + hi + 
h r = h v>r — htp t i = AIt, where A is constant in time and 
It is a tensor product of the identity operators of the indi- 
vidual islands. The conservation reduces to a comparison 
of the shunt capacitors since the junction charges cancel. 
This means that if some charge is added to one shunt 
capacitor, an equivalent charge must be removed from 
the other. Hence, we can rewrite + Al T 

and n = = - h Vl i - 4 I T - We thus ob- 

tain the CPN operators for the left and right islands as 
m = h - | + A\ T and n r = — n — § + ^It- Note that 
for the charge transport, the total charge in the system 
is irrelevant and can be taken to zero without loss of 
generality, that is, h v _ r — n Vl i = AIt = 0. We adopt 
this convention. In addition, we omit explicitly writing 
all identity operators in the following presentation of the 
Hamiltonian. 

Deriving the charging Hamiltonian for the device is 
a straigthforward task using the CPN basis formulation 
for the system of islands. After applying the above- 
mentioned transformations, the charging Hamiltonian as- 
sumes the form 

H ch = E c (h-n g ) 2 + E v h 2 + H FIN , (Al) 

where we have 



E, 



2e' 



2ef 



(2C V + C g )(2C v + Ci + C r ) 



^{C g ClC r + CUC g +Q + C r ) + C v [2ClC r + C g (Cl + C r )}} ' 



C g (Cl + C r ) + 2C v (C g +Cl+ C r ) 



(A2) 



CgClCr + CUC g + Cl+ C r ) + C v [2ClC r + Cg{Q + C r )} 
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Here C g is the gate capacitance, C v is the capacitance final term in the Hamiltonian is 
of one of the identical shunt capacitors and C\ , C r are 

the internal capacitances of the tunable junctions. Fur- #fin = E\n g h + E 2 n 2 + E 3 n g h + E^nn, (A3) 

thermore, the normalized gate charge is n g = V g C g /(2e), 

where V g is the externally controlled gate voltage. The where the prefactors are 





= 2e 2 


E 2 


= 2e 2 


E, 


= 2e 2 


Ei 


= 2e 2 



C g (2C ip + Ci + C r ) 



2{C g dC r + C 2 (C g +Q + C r ) + C v [2dC r + C g (Ci + C r )}} : 

4QCr - CgjCl + C r ) + 2C v {d + C r - C g ) 
4{C g C t C r + C 2 JC g +Q + C r ) + C v [2dC r + C g (d + Cr)}}' 



2C ip (C r — C{) 



(A4) 



C g ClCr + C 2 (C g +d+ C r ) + C v \2dC r + C g (d + C r )] ' 

(2d v + C g )(d - C r ) 

CgClCr + C 2 (C g +d + C r ) + C v \2dCr + C g (d + C r )} ' 



Presenting the Hamiltonian in this form showcases 
its properties. If the island is capacitively symmetric 
C; = C r = Cj, we have £3 = E4 = so that the 
first-order feed term and the term coupling the island 
and feed degrees of freedoms vanish from the charging 
Hamiltonian. Furthermore, the remaining prefactors be- 
come En - 2e 2 (^+c g )(c v + Cj ) 

come G C 2 [C g C*+Cl(C g +2Cj)+2C v Cj{Cj+C a )V 
2 2 2(C g Cj+C v (C a +2C.,)) 



Ey 

Ei 



E> 



C g C 2 J +C2(C g +2Cj)+2C v Cj(Cj+C g ) ' 



2e 2 



C g Cj+Cl{C g +2Cj)+2C v Cj(Cj+C g ) 



2e 2 



2Cj-C g Cj+C^(2Cj-C g ) 



and 
If we 



2[C g C 2 J +Cl(C g +2Cj)+2C f Cj{Cj+C g )Y 

assume that the gate capacitance is vanishingly small 
compared to the shunt capacitance, the first order island 
charge term becomes negligible. The remaining term 
E 2 n 2 in J?fin is proportional to the identity operator, 
and hence only induces a time-dependent shift of the 
zero point of the energy that can be neglected. 

The inductive part of the Hamiltonian stems from the 
inductive energy Hi n d = \EI\\ where L is the total in- 
ductance of the loop element and II is the operator for 
the current passing through it. Note that in addition 
to the geometric inductance Lq 1 L = Lq + Lk also 
includes the kinetic inductance Lk which can poten- 
tially be significant for a superconducting wire^. The 
total magnetic flux threading the loop can be expressed 
as a sum of the applied external flux and the flux in- 
duced by the supercurrent due to the geometric induc- 
tance $ = <f> C xt + LqIl- Furthermore, the operator 
for the total phase difference over all the weak links 
in the loop (p is related to the total flux operator by 



^■L K I L =<p- ^-L K I L , where (pi and 
fir are the operators for the gauge-invariant phase differ- 
ences over the left and right SQUIDs, respectively. The 
kinetic inductance adds to the total phase difference over 
the loop according to the ac Josephson equation^, and 
the directions over which the phase differences and the in- 



ductive current are calculated are given in Figs. [T] and [5] 
Using the above-mentioned identities, the operator for 
the inductive current assumes the form 



f $0 v 



¥0), 



(A5) 



where ipo = 2?r ( |' cxt and we denote the flux quantum by 
$0 = h/(2e). Thus, the inductive Hamiltonian can be 
written as 



#ind = -£-{$- ^o) 2 , 



(A6) 



L\2ty) 



where we denote El 

Finally, the Josephson part describing the tunneling 
through the junctions is given by^ 



Hj = -E Jr COS (p r - Eji COS Ifii, 



(A7) 



where Eji($i) and Ej r ($> r ) are the Josephson energies 
for the left and right tunable junctions, respectively. 
Starting from the operators for the phase differences 
over the SQUIDs, we perform a transformation to a new 
pair of phase operators such that <p — (pi + tp r and 
4> = (if r — This transformation is well-founded 

since the canonical commutation relation between the 
charge and phase operators of the individual junctions 
yields [4>, h] = i indicating that <f> is the island phase 
operator and [ip,n] = i indicating that ip is the canoni- 
cal conjugate of the feed CPN operator—. If the average 
charge imbalance between the left and right leads changes 
by one, the occupation in the eigenbasis of n subsequently 
changes. This implies charge transfer in the system and 
allows us to monitor the system currents. With the trans- 
formed phases, the Josephson part assumes the form 



I - 



(A8) 
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Appendix B: Perturbative analysis H$ and e k is the corresponding eigenenergy. Denoting 

the adiabatic states and energies as |$fc rl ) and Ek n corre- 
The unperturbed states \k) |* n ) are defined via the sponding to the fcnth unperturbed state, the inner prod- 



ei genprob lem [H s ® I v + I„ (8 H v ] \k) |*„) = [e k + 
JlE v Ei, (n + hj] \k) \^ n ) where \k) is an eigenstate of 



ucts relating to the different current elements mentioned 
in Sec. MI Al become 



($ fcn |/j 0) |$ pm ) 



{k\I S \p)S n , m + go-^2(k\I s \l) (1\I S \P) 



— £/ + y / 2E tp EL — — y / 2E ip Ei J 



m — 1 



VmTT 



+ 



e p — £/ — \/2E v El efe — e/ + ^J2E v El 



m+1 



O(.9o 2 ), 



(Bl) 



and 



($kn\Ii 1] \$pm) = -ga^r (k\[H B + Hl]\p) (y/mS^m-i + y/m + l<J„, m+ i) + 0{gl), 



(B2) 



and 



<f>kn\h\<5> pm ) = (k\I s \p) 



2E v El 



(e P ~€ k f- 2E V E L 



S n , m + j(%E v El)*5 ktP {^i8 n , m -i + Vm+l<J„, m +i) + 0(sg). (B3) 



Using Eqs. (|B1|) and (|B2[) . the dynamical current through Eq. (|B3[) allows us to write I L ' n> = ($fc„|Iz,|<I>fen) = 
the island carried by |$fc„) can be readily written as _ (k\i s \k) + 0(<7q). 

j(D,kn) _ ^$^1^1$^^ = (k\Is\k) + O(g^). Similarly, In order to study the geometric currents, we need the 

temporal derivatives of the adiabatic states. The corre- 
sponding inner product is 



($pm\d t \®kn) = (P\k) d n . rn + g - ^ 



(l\k) Ws\l) 



go-ot 



w) m\k) 

m5n,m-i 



mSr, 



— e; + \/2E v El — e/ — ^J2E v El 
y/m + 16 nm+ i 



k £ p — e z + \/2E v El e p — £; — yj2E tp Ei J 
Ws\k) 



(B4) 



Vm + l<y„, m+ i +go-d t - 
e \e k - t P + \J2E ip E L e V £fe - e p - ^2E V E L 



md n ,m-l 



O(g 2 o), 



where the dot indicates time derivative. Using Eq. ()B4|) 
and the above-mentioned current elements, the geometric 
current through the island becomes 
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^kn ^pm , , ^A: 



where ^ denotes a double summation that does not 
include the term for which p = k and m = n simulta- 
neously. Calculating the inductive element in a similar 
manner, we obtain jf ,kn ' = _/j G > fc ") U p ^ G ^} ie g rs ^ 
order. 



Appendix C: Parameter cycle 

The parameter cycle consists of six equivalently long 
legs designed to maintain good charge quantization per 
cycle. In the beginning, both SQUIDs are closed so that 
their residual Josephson energies are E") m and i?""", and 
the gate charge is n™ m . Nonzero values for the residual 



Josephson energies are typically used to include the effect 
of nonidealities stemming from the manufacturing of the 
junctions, and lead to a nonvanishing supercurrent flow- 
ing through the device. During the first stage, the left 
SQUID is opened by bringing its Josephson energy lin- 
early to Ejf*. During the second stage, the gate charge 
is linerly increased to n™ ax drawing charge from the left 
lead to the centre island. The left SQUID is then linearly 
closed and the right SQUID opened in a similar manner 
bringing its Josephson energy to Eff*. During the fifth 
stage, the gate charge is brought linearly back to n™ ln 
pushing the charge to the right lead. Finally, the right 
SQUID is linearly closed and the system has completed 
a closed loop in the external parameter space. 
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